Massive phonon modes from a BEC-based analog model 
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Two-component BECs subject to laser-induced coupling exhibit a complicated spectrum of exci- 
tations, which can be viewed as two interacting phonon modes. We study the conditions required 
to make these two phonon modes decouple. Once decoupled, the phonons not only can be arranged 
travel at different speeds, but one of the mod es can be g iven a mass — it exhibits the dispersion 
relation of a massive relativistic particle u) = yWo + c 2 k 2 . This is a new and unexpected excitation 
mode for the coupled BEC system. Apart from its intrinsic interest to the BEC community, this 
observation is also of interest for the "analogue gravity" programme, as it opens the possibility for 
using BECs to simulate massive relativistic particles in an effective "acoustic geometry" . 
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Analogue models for gravitation can be used to simu- 
late classical and quantum field theory in curved space- 
time. The first analogue model for black holes, and for 
simulating Hawking evaporation was suggested by Bill 
Unruh [fj. He demonstrated that a sound wave propa- 
gating though a converging fluid flow exhibits the same 
kinematics as light does in the presence of a curved space- 
time background. Since then several other media have 
been analyzed, and the field has developed tremendously. 
The first approach specifically using Bose-Einstein con- 
densates as an analogue model was made nineteen years 
later . Since then various configurations of BECs have 
been studied to simulate different scenarios for grav- 
ity EJ El H- Until now it has only been pos- 
sible to simulate light, (generally speaking, massless rel- 
ativistic particles) , propagating through a curved space- 
time 0, 0, 0, . In the following a two-species 
BEC is used to extend the class of equations that can be 
simulated to the full curved-space Klein-Gordon equa- 
tions. In the language of the BEC community — we 
have developed a way of giving a mass to the phonon. 

The system we will use in our theoretical analysis is an 



ultra cold 2-component BEC atomic gas. For example, a 
two-component condensate of 87 Rb atoms in different hy- 
pcrfmc levels, which we label \A) and \B). (Experiments 
using two different spin states, \F = 1, m = — 1) and If = 
2,m = 2), were first performed at JILA in 1999 14].) 
At zero temperature nearly all atoms occupy the ground 
state. Then the quantized field describing the micro- 
scopic system can be replaced by a classical mean- field, 
a macroscopic wave-function. In this so-called mean- 
field approximation the number of non-condensed atoms 
is small. Interactions between the condensed and non- 
condensed atoms are neglected in the mathematical de- 
scription, but two-particle collisions between condensed 
atoms are included. In the case of a two-component sys- 
tem, interactions within each species (Uaa, Ubb) and 
between the different species (Uab = Uba) take place. 
In addition the two condensates are coupled by a laser- 
beam, which drives transitions between the two hyperfine 
states with a constant rate A. (Without the coupling A 
no mass term is generated, which is consistent with |8(.) 
The resulting coupled time-dependent Gross-Pitaevskii 
equations are: 



ihd t ^A 
ihdt^B 



V 2 + VX + Uaa \*a\ 2 + U AB |*s| 2 



2 rriA 
Ivhb 



V 2 + Vi + U BB \^b\ 2 + U AB \^a\ 2 



f A + A* b , 



(1) 



where VXjB = ^a/b ~ I^A/b denotes the combined effects 
of the external potential V A IB an d the chemical poten- 
tial Ha/b |l5lll6|. In the eikonal approximation the two 
stationary background states are described by their den- 



sities {pa,Pb} and phases {Oa, 9b}' 

= ^/px e l{6x] for X = A, B . (2) 

These four variables are in general not independent of 
time and space. 
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In the following, we study zero sound in the overlap 
region of the two-component system, produced by excit- 
ing density perturbations which are small compared to 
the density of each condensate cloud. In the first exper- 
iment studying localized excitations in a one-component 
Bose-Einstein condensate [l^, a laser beam was used 
to generate a small modulation in the density. Using a 
phase-contrast imaging it was shown that the resulting 
perturbation corresponds to a sound wave. The observed 
speed of sound is 



c(r) 



4-7T H 2 a p(r) 



Up(r) 



(3) 



where p(r) is the density of the ground state, a is the 
scattering length, m is the atomic mass and U is the 
self-interaction constant. The mathematical equations 
describing these perturbation leads to the well-known 
hydrodynamic equations, which are the basis for the 
most fruitful of the analogies between condensed matter 
physics and general relativity Q, 0, 0, • 

The same method can be used to obtain the kinematic 
equations for small perturbations propagating in a two- 
component system. Given that the density modulation 
is small, the perturbations in the densities and phases 
can be linearized around their unperturbed macroscopic 
states {pao,Pbq}- 

*x = Vpxo + epxie l{9xo+Eexl) for X = A, B . (4) 



These states still satisfy the coupled Gross-Pitaevskii 
equation. After a straightforward calculation the terms 
of first order in e representing the sound waves include 
two coupled equations for the perturbation of the phases 



_ V7Q U A a . U A b 

VAOVVAI r— PA1 H r— PB1, 



'Bl 



vbo^Obi - — |^ PBl 



h 

Uab 



(5) 



' PAl- 



Here vao = Oi/ttia^Oao and vbo = (fi/ 1713)^ 8 bo are 
the background velocities for the macroscopic ground 
states, while 



Uaa 
Ubb 
Uab 




(0) 



2 ^/Pao Pbo ' 



are modified interaction potentials for the two coupled 
condensates. In addition to these two phase equations, 
there are two coupled equations for the density pertur- 
bations 



, h ^ \ 2A , . 

PAl = - V ( PAO^VAl + PAlVAO I + -ry/PAO PBO {?B1 - &AX), 

k t<ia / n 

Pbi = - V ( pbo^Obi + pbivbo J + ~r \Jpaq Pbo (Oai - Obi)- 

1 n%B In 



(7) 
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It is useful to define the coupling matrix 

E = -( ^ AA ^~ AB 
H V -Uab U B b 

A second coupling matrix can be introduced as 



A = 



2y/PA0 PBO ( +A -A 

-A +A 



(8) 



(9) 



which vanishes completely if the coupling laser is 
switched off. Last but not least, it is also useful to intro- 
duce the background velocity matrix V (a 2 x 2 matrix 
of 3- vectors) 



V = 



VAO 
VBO 



(10) 



and the mass-density matrix D 
D = H 



PAO Q 

Q Pbo 
m B 



(11) 



Collecting terms into a 2 x 2 matrix equation, the equa- 
tions for the phases JSJ and densities Q become 



V ■ V6, 



(12) 



p= -V (D- S79 + V ■ p) - A- 6, (13) 

where 9 T = {6ai,0bi) and p T = (pauPbi)- 

Equation l|12|) can be used to eliminate p and p in 
equation (|13l) . leaving us with a single matrix equation 
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for the perturbed phases: 

dt (V 1 • 0) = -a t (s- 1 ■ V ■ ve) - v (v ■ s- 1 • 

+V [(£> - V • S" 1 • V) V0] + A-0 (14) 

This equation tells us how a localized collective excita- 
tion in a two-component system with a permanent cou- 
pling — if A is nonzero — develops in time. So far there 
are no restrictions on the masses (m^ms), densities 
(pao, Abo)) background velocities (vaojVbo), interaction 
constants (Uaa,Ubb,Uab), and coupling constant (A). 
In the eikonal approximation this differential equation 
leads to the Fresnel equation 

det {lo 2 S- 1 - w [S- 1 • (Vk) + (Vk) • S- 1 ] 

-[D k 2 -(Vk)-E- 1 -(Vk)] + A}}= 0,(15) 

which is in general a quartic dispersion relation for two 
interacting phonon modes. Here Vk denotes the matrix 



Vk = 



VAo-k 
vbo ■ k 



(16) 



The first step in analyzing equation (|14J) is to ask 
whether it is possible to decouple the system into two 
independent phonon modes. We have found decoupling 
is not possible without introducing several constraints on 
the background quantities. 

Focusing on the last term in equation (|14(1 . the eigen- 
vectors for a non-zero coupling A 7^ are given by 
{[1, 1], [— 1, 1]} and the corresponding eigenvalues are 
{0, 4 XytpAoPBo/f 1 }- The eigenvectors are fixed and in- 
dependent of any of the other physical variables. As a 
result the only way to decouple equation (|14|l into two 
phonon modes, is to decompose it in the following way: 



(17) 



We now analyze equation l|14|) term by term with respect 
to this decomposition. 

The term on the LHS has the same eigenvectors as 
equation ((T7|l if and only if Uaa — Ubb- The eigenvalues 
of corresponding to {[1, 1], [—1, 1]} are {ft/ (Uaa — 
Uab), ^/(Uaa + Uab)}- This places another constraint 
on the interaction variables: Uaa 7^ Uab- (Further dis- 
cussion regarding this tight requirement will be presented 
later.) 

The eigenvectors for the first two terms on the RHS in 
equation (jl4|) can be found in one step, by using the fact 
that they have simultaneous eigenvectors if and only if 
the commutator [S _1 , V] — vanishes. Therefore vao — 
vbo = Vq, the background velocities must be equal if the 
two phonon modes are to decouple. The corresponding 
eigenvalues are those from the matrix multiplied by 
Wo- In other words, the backgrounds of two condensates 
must be in phase 0ao = Obo- 



We are now left with the penultimate term in equation 
(fHfjl . Because V = vq ■ 1, and the eigenvectors of S _1 are 
already known, there is only the mass-density matrix D 
to consider. The last constraint to decouple the equation 
for the two phases then is HpAo/ m A — frpBo/mB = d. 
So D = d ■ 1 has for both eigenvectors the eigenvalue 
d = hy/pAo PBo/^m A m B - 

Applying all this to equation (fl4*jl one obtains two de- 
coupled equations for the phonon modes described by the 
eigenstates (|T7j) : 



621 X , for 7=1,2. 



Here 



(18) 



(19) 



where the propagation speeds are defined in terms of the 
eigenvalues Sj of the matrix H by 



=? = Sr d 



d{U A A + (-iyUAB) 

h 



(20) 



Introducing effective "spacetime metrics" by the iden- 
tifications ^J—gi gj V — fj" and gi = 1/ det^j"], we can 
recast these equations as a pair of curved-space Klein- 
Gordon [massive d'Alembertian] equations 



gigfa v 0i)-A J 6 3I I = O. (21) 



Here 



9i 



-2/(23-2) 





~ v 


-v> 


Cj5 13 — Vq Vq 



which depends on the space-time dimension D 
such a manner that 



C = 1 - 



22) 
in 



2) 




-V 3 Q 









Finally the mass-term is 



^ = J_WPMPm (c 2 /d 23 } V(Z3-2) 



corresponding to the natural oscillation frequency 
2/(23-2) , > , 2 



(23) 



(24) 



w„ = m 2 Cn 



hd 



(25) 



The dispersion relation (in the eikonal limit) is then 

(w - v„ ■ k) 2 - cjk 2 = lo 2 S 2I . (26) 
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(A similar calculation, but restricted to a one-condensate 
system, where all variables are allowed to be time and 
space dependent, but no mass term is present, has been 
presented in 0.) Equation (|2T)l is the curved-space 
Klein-Gordon equation (massive d'Alembertian equa- 
tion). For 6 1, (corresponding to perturbations in the two 
condensates A and B oscillating "in phase"), the mass 
term is always zero. However, for a laser-coupled system 
(A 7^ 0) the mass-term in the equation for 62, (corre- 
sponding to perturbations in the two condensates A and 
B oscillating "in anti-phase"), does not vanish. 

Comparing the definition for the speed of sound in 
a one component system, with the cj introduced here, we 
see that the cj ifflty are the modified speeds of sound for 
each phonon mode. (If U ab = A = the two condensates 
decouple and the c/ limit to the independent phonon 
speeds in each condensate cloud.) This fact leaves us 
with the possibility of constructing two different types of 
analog model. So far we have been dealing with a two- 
metric structure, which is interesting in itself 0, llij . 
For instance, in the absence of laser-coupling (A = 0) the 
presence of two different speed of sounds can be used for 
tuning effects 0- 

If we wish to more accurately simulate the curved 
spacetime of our own universe, another constraint should 
be placed on the system, to make the two speeds of sound 
equal c\ — ci- This yields a single sound-cone structure, 
to match the observed fact that our universe exhibits a 
single light-cone structure. This condition is fulfilled if 
we set U ab — 0. While the in-phase perturbations will 
propagate exactly at the speed of sound, 

v s = v + k CI , (27) 

the anti-phase perturbations will move with a lower 
group velocity given by: 

3k ^ § 2I + c 2 k 2 

Here k is the usual wave number. This explicitly demon- 
strates that the group velocity of the anti-phase eigen- 
state depends on the laser-induced coupling between the 
condensates. 

In conclusion, the calculation presented in this article 
is of interest to two separate communities. For the BEC 
community, it provides a specific example of how to tune 
an interacting 2-BEC condensate in such a manner as 
to obtain a massive phonon. Without the fine tuning, 
it provides an example of two interacting phonon modes 
whose dispersion relation is governed by the quartic Fres- 
nel equation I] 15(1. For the general relativity community, 
this article provides an example of an analogue system 
that can be used to mimic a minimally coupled scalar 
field embedded in a curved spacetime. 
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